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Abstract — Temporal and energy characteristics of the very-high-energy gamma-ray bursts from evaporating primordial 
black holes have been calculated by assuming that the photospheric and chromospheric effects are negligible. The 
technique of searching for such bursts on shower arrays is described. We show that the burst time profile and the array 
dead time should be taken into account to interpret experimental data. Based on data from the Andyrchy array of the 
Baksan Neutrino Observatory (Institute for Nuclear Research, Russian Academy of Sciences), we have obtained an upper 
limit on the number density of evaporating primordial black holes in a local region of space with a scale size of ~ 10 -3 
pc. Comparison with the results of previous experiments is made. 
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Introduction 

Primordial black holes (PBHs) can be formed in the 
early Universe through the gravitational collapse of 
primordial cosmological density fluctuations. There- 
fore, the formation probability of PBHs and their ob- 
servational manifestations depend significantly on how 
the primordial density fluctuations emerged and devel- 
oped. The pattern of black hole formation is deter- 
mined not only by the cosmology and physics of the 
early Universe. Theoretical predictions of the PBH 
formation probability depend strongly on the adopted 
theory of gravitation and, which is also important, 
on the adopted model of gravitational collapse. The 
evaporation of black holes on which their experimen- 
tal search is based has not been completely studied 
either. Thus, PBH detection will provide valuable in- 
formation about the early Universe and can be a unique 
test of the general theory of relativity, cosmology, and 
quantum gravity (Carr 2003). Knowledge of the spa- 
tial distribution of PBHs is important for their direct 



search. As was shown by Chisholm (2006), the local 
PBH number density in our Galaxy could be many or- 
ders of magnitude higher than the mean PBH number 
density in the Universe (the density ratio could reach 
~ 10 22 , which is much larger than the previously pre- 
dicted value of ~ 10 7 ). Therefore, the constraints on 
the PBH number density imposed by direct searches 
can be more stringent than those imposed by diffuse 
extragalactic gamma-ray background measurements. 

The high-energy gamma-ray bursts (GRBs), i.e., 
significant and time-localized excesses of gamma radi- 
ation above the background, are generated at the final 
PBH evaporation stage. Since the calculated tempo- 
ral and energy characteristics of such bursts depend on 
the theoretical evaporation model (Bugaev et al. 2007) , 
the technique of an experimental search and, accord- 
ingly, the constraints imposed on the PBH number den- 
sity in the local Universe are modeldependcnt. PBHs 
can be searched for in experiments on arrays designed 
to detect extensive air showers (EASs) from cosmic 
rays with effective primary gamma-ray energies of 10 
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Figure 1: Probability of EAS detection by the 
Andyrchy array versus primary gamma-ray photon en- 
ergy for various zenith angles 9. 

TeV or higher only within the framework of the evapo- 
ration model without a chromosphere (MacGibbon and 
Webber 1990). The duration of the high-energy GRBs 
predicted by chromosphcric evaporation models is too 
short, much shorter than the dead time of EAS arrays. 

It should be noted that the duration of the high- 
energy GRBs is fairly short in the evaporation model 
without a chromosphere as well. Therefore, the effect 
of the array dead time on the burst detection proba- 
bility should be taken into account when interpreting 
the experimental data from EAS arrays with a high 
threshold energy of the primary gamma-ray photons. 

The Experiment 

The Andyrchy array to detect EASs from cosmic rays is 
located on the flank of Mount Andyrchy at an altitude 
of <~ 2060 m above sea level; its geographical coordi- 
nates are 43.28° N and 42.69° E. The array consists of 
37 scintillation detectors based on plastic scintillator, 
each with an area of 1 m 2 . The separation between 
the detectors in the horizontal plane is ~ 40 m and 
the total area of the array is 5 x 10 4 m 2 . The EAS 
trigger becomes active when > 4 array detectors are 
triggered simultaneously; the trigger rate is ~ 9 s _1 . 
The effective angular resolution of the array for such 
events is 3.8°. The array dead time per EAS event does 
not depend on the EAS power and is 1 ms. The array 
and its operating parameters were described in detail 
previously (Petkov et al. 2006). 

The detection probabilities P{E~ n 8) of the EASs 
generated by primary gamma-ray photons with energy 
E 1 incident on the array at zenith angle 9 were deter- 
mined by simulating electromagnetic cascades in the 
atmosphere and the detector (Smirnov et al. 2005; 
Smirnov 2005). The CORSIKA code (version 6.03, the 
QGSJET model) (Heck et al. 1998) was used to simu- 



late electromagnetic cascades in the atmosphere. The 
characteristics of the secondary particles that reached 
the array level were used as input parameters in the 
code for calculating the detector response, in which the 
energy release and the triggering time of each array de- 
tector were calculated. For the event simulated in this 
way, the arrival direction of the simulated EAS was re- 
constructed on the basis of a standard technique used 
in processing the recorded events. In Fig. 1, the EAS 
detection probability is plotted against the primary 
gamma-ray photon energy for several zenith angles. 
Since the detection probability of primary gamma-ray 
photons is a relatively smooth function of the photon 
energy, the median energy of the primary gamma-ray 
photons detected by the array depends on their energy 
spectrum. Following Alexandreas et al. (1993), we 
will take the median energy of the primary gammaray 
photons when the source is located at zenith and the 
gamma-ray spectrumis a power law with an index of 
—2.7 as the effective energy of the gamma-ray photons 
detected by the array. For the EASs detected by the 
Andyrchy array, this energy is 60 TeV. 

For each of the events selected by the EAS trig- 
ger, we reconstructed the EAS arrival direction, i.e., 
the zenith and azimuth angles (9, <j)) in the array co- 
ordinate system. Based on our processing, we created 
an archive of preprocessed information for the period 
1996-2001 (the net accumulation time is ~ 1100 days 
and the total number of events is ~ 6.22 x 10 8 ), which 
contains the absolute event time (with an accuracy of 1 
ms) and the EAS arrival direction. Searching for GRBs 
over the celestial sphere (without referencing to the al- 
ready detected bursts) is, in fact, searching for spa- 
tiotemporal concentrations of events (clusters). Since 
we take fairly short time intervals, spatial concentra- 
tions of events are searched for in the horizontal coor- 
dinate system. For each event i with an absolute time 
ti and arrival angles (9,<p)i, we search for a cluster of 
such events i , i + 1 , . . . , i + n — 1 that the EAS arrival di- 
rections differ by less than a r from the weighted mean 
direction. Thus, each cluster is characterized by mul- 
tiplicity n, duration At, absolute time T, and arrival 
direction [9, <j>). 

Previously, data from the Andyrchy array (Smirnov 
et al. 2005; Smirnov 2005) were used to search for cos- 
mic GRBs over the celestial sphere. Groups of EASs 
arrived from one angular cell a r = 4.0° in radius were 
selected; the minimum and maximum time differences 
in the cluster were taken to be 10 ms and 10 s, re- 
spectively. For each multiplicity (n > 2), the depen- 
dences of the number of clusters with a given multiplic- 
ity on the interval At were derived. The background 
of chance coincidences (the formation of clusters with 
a given multiplicity n) was calculated using a similar 
processing of the simulated events. The EAS arrival 
angles (9, <p) and the time between the EAS arrivals 
were simulated using experimental distributions. The 
experimentally measured detection rates of such clus- 
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ters are in agreement with those expected from the 
background of chance coincidences (within the limits 
of one standard deviation) . 

Searching for GRBs from PBHs 
for the Evaporation Model With- 
out a Chromosphere 

The spectrum of the gamma-ray photons emitted by 
PBHs, dN^/dE 7 , depends on the time ti left until 
the end of black hole evaporation. In the evapora- 
tion model without a chromosphere (MacGibbon and 
Webber 1990), the numerically calculated and time- 
integrated photon spectrum can be fitted by a piece- 
wise power function (for a black hole temperature 
Th 3> and gamma-ray energies E 1 > m^/2): 



3/2 



dN 7 

~dK 



(*).(*) 



£-7 < Th 



(Bugaev et al. 2007), where Eq = 10 5 and all energies 
are measured in GeV. The parameter N is 



N rj VQ x 7- 10 19 , 

while vq includes the effective number of degrees of free- 
dom of the quarks and gluons used in calculating the 
photon spectrum. A simple estimate for N includes 
the u and d quarks, their antiquarks, and all gluons; 
in this case, v = (1/3) (3 • 2 • 2 • 2 + 8 • 2) = 40/3 (1/3 
is the ratio of the number of ir° mesons to the total 
number of ir mesons). The corresponding value of No 
is 9 x 10 20 . The lower limit for the parameter v$ can 
be obtained by taking into account only one type of 
quarks {y^ = 4); in this case, Nq w 3 x 10 20 . 

The time until the end of PBH evaporation ti (in 
seconds) is related to the black hole temperature Th 
(in GeV) by 

ti 



4.7 x 10 11 T H 3 . 



Since the array detects the EASs generated by pri- 
mary gamma-ray photons with energy E 1 incident at 
zenith angle 9 with probability P(E 7 , 9) and since the 
spectrum of the gamma-ray photons emitted by PBHs 
dN^/dEy depends on ti, the spectrum recorded by the 
array (its response function) P(E 7 , 9)dN 1 / dE 1 also de- 
pends on ti. Here, 9 is the zenith angle at which the 
PBH is seen from the array. Figure 2 shows the photon 
spectra and response functions of the Andyrchy array 
for the vertical direction and for ti — 1 s and 1 ms, 
with N = 3 x 10 20 . 

The total number of gamma-ray photons emitted 
by PBHs that can be detected by the array, 
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Figure 2: (a) Photon spectrum from evaporating 
PBHs. (b) Response functions of the Andyrchy array 
for the vertical direction at two values of U. 



also depends on t\. Figure 3 shows the dependences 
N 7 (6, U) (integrated burst profile) for several zenith an- 
gles 9. It should be noted that the burst profile mainly 
determines the technique of searching for GRBs from 
PBHs on a specific array, in particular, the choice of a 
time interval At for the search for EAS clusters. Let 
us call the time left until the end of black hole evap- 
oration during which 99% of the gamma-ray photons 
that can be detected by a given array are emitted the 
burst duration % for this array. In Fig. 4, t\> is plot- 
ted against the zenith angle for the Andyrchy array. 
We see from the figure that the array can record the 
events from evaporating PBHs in a limited range of 
zenith angles, because the expected GRB duration at 
large zenith angles is shorter than the array dead time 
per event td- 

We searched for the groups of EASs arrived from 
one angular cell for time intervals of a given duration, 
At = tb(0). Because of the short time intervals (accord- 
ingly, the reduced background of chance coincidences), 
we used large angular cells: a r = 7.0°; such a cell con- 
tains 90% of the events from a point source. The dis- 
tributions of the detected clusters in multiplicity were 
obtained for At = 40, 24, 11, and 3 ms, which corre- 
spond to the zenith angles 6 = 0°,10°,20°, and 30°. 
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Figure 3: Total number of gamma-ray photons emit- 
ted by PBH that can be detected by the Andyrchy 
array versus the time ti left until the end of black hole 
evaporation: 1, = 0°; 2, = 10°; 3, = 20°; 4, 
= 30°. 



Figure 4: Burst duration versus zenith angle for the 
Andyrchy array. 



and the probability to detect all n EASs over the burst 
duration tb with the array dead time per event td- 



The results of our search for two values of t b are pre- 
sented in Fig. 5. The maximum multiplicities n'(0) of 
the detected clusters are 4, 4, 3, and 2, respectively. 

Results 

We obtained a constraint on the number density of 
evaporating PBHs using the technique described in 
Linton et al. (2006). Let a PBH be located at distance 
r from the array and be seen from it at zenith angle 0. 
The mean number of gamma-ray photons detected by 
the array over the burst duration tb is then 



F(n,n,t b ) = f(n,t b ,t d )- 



n{0) = 



eN(t b (0))S(0) 
Anr 2 



where S(0) is the array area and e — 0.9 is the frac- 
tion of the events from a point source that fell into an 
angular cell with a size of 7.0°. The number of bursts 
detected over the total observation time T can be rep- 
resented as 

N = PpbhTVeff , 

where 

oc 

V e/J = [dSlf drr 2 F(n(0),n(0),t b (0)) 



is the effective volume of the space surveyed by the 
array, p pb h is the number density of evaporating PBHs, 
and F(n, n, tb) is the detection probability of a cluster 
of n EASs with the mean value of n over the burst 
duration %■ In turn, the total detection probability can 
be expressed as a product of the Poisson probability of 
n EASs falling on the array with the mean value of n 



In calculating the effective volume, we take n{0) = 
n'(0) + 1 (this means that the distributions of the de- 
tected clusters in multiplicity can be explained by the 
background of chance coincidences). The probabilities 
f(n, tb, td) of detecting n{0) EASs over the burst dura- 
tion tb{0) with the array dead time per event td were 
calculated by the Monte Carlo method by taking into 
account the burst time profile for a given zenith angle. 
The effective volume of the space surveyed by the ar- 
ray, V eff , is 1.88 x 10~ 9 pc 3 , with N Q being 9 x 10 20 . If 
the evaporating PBHs arc distributed uniformly in the 
local region of the Galaxy, then the upper limit pi im on 
the number density of evaporating PBHs at the 99% 
confidence level can be calculated from the formula 



Pi i 



4.6 



VeffT ' 



m our case, pu r , 



8.2 x 10 8 pc~ 3 yr" 



Comparison with the Results of 
Previous Studies and Conclusions 

Previously, the high-energy GRBs from evaporat- 
ing PBHs were searched for in experiments on the 
CYGNUS (Alexandreas et al. 1993; the threshold en- 
ergy of the primary gamma-ray photons is E t h ~ 30 
TeV), HEGRA (Funk et al. 1995; E th - 30 TeV), 
and Tibet (Amcnomori et al. 1995; E th - 10 TeV) 
EAS arrays for fixed one-second time intervals with- 
out taking into account the burst duration (and its 
dependence on the zenith angle for the specific array). 
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Figure 5: Measured number of clusters versus cluster 
size for two time intervals, At — 40 ms and At = 11 
ms. 

In these experiments, constraints on the number den- 
sity of evaporating PBHs were obtained by assuming 
the absence of any fluctuations in the number of EASs 
falling on the array (it was assumed that n — n) and 
by disregarding the effects of the burst time profile and 
the array dead time on the burst detection probability. 
The latter is particularly important for arrays with a 
high gamma-ray detection threshold (i.e., with a short 
burst duration U) and a relatively long dead time td- 
Disregarding the effects of the array dead time and the 
burst time profile on the burst detection probability 
when interpreting the data will lead to an appreciable 
overestimation of this probability (i.e., the factor F in 
the above expression for the effective volume V e ff) and, 
as a result, to an underestimation of the upper limit 
on the PBH number density. For example, in this case, 



a constraint of 3.1 x 10 6 pc -3 yr -1 instead of the real 
8.2 x 10 8 pc~ 3 yr _1 is obtained for the Andyrchy ar- 
ray. In the limiting case, the array can be insensitive 
to the high-energy GRBs from PBHs altogether, since 
for td > ti, no more than one gamma-ray photon can 
be detected over the burst duration and a burst of n 
gamma-ray photons can be detected if ntd < 

Figure 6a presents the results of the studies listed 
above as they are given by their authors. This fig- 
ure also shows the constraint obtained at ~ 1 TcV 
in an experiment on the Whipple Cherenkov telescope 
(Linton et al. 2006) for a time interval of 5 s. In addi- 
tion, for comparison, it also shows our constraint that 
we obtained in the approximation of a zero array dead 
time and by disregarding the burst time profile. The 
effective energy of the detected gamma-ray photons is 
along the horizontal axis. Two constraints obtained 
for different event selection conditions are given for the 
HEGRA array. 

The upper limits on the number density of evapo- 
rating PBHs with allowance made for the burst time 
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Figure 6: Upper limits on the number density of evap- 
orating PBHs for the evaporation model without a 
chromosphere versus effective energy of the detected 
gamma-ray photons. Panel (a): 1, the Andyrchy ar- 
ray (this paper, without the effects of the burst time 
profile and the array dead time on the burst detec- 
tion probability); 2, original data from the Tibet array 
(Amenomori et al. 1995); 3, the Whipple Cherenkov 
telescope (Linton et al. 2006); 4, original data from the 
CYGNUS array (Alexandreas et al. 1993); 5, original 
data from the HEGRA array (Funk et al. 1995). Panel 
(b): 1, the Andyrchy array (this paper, with the effects 
of the burst time profile and the array dead time on 
the burst detection probability); 2, recalculation of the 
Tibet data (see the text); 3, the Whipple Cherenkov 
telescope (Linton et al. 2006). 



profile and the array dead time are shown in Fig. 
6b. We obtained the constraint for the Tibet array 
(7.6 x 10 6 pc~ 3 yr -1 ) by recalculating the data from 
Amenomori et al. (1995) using the technique described 
above (the dead time of this array is 5 ms per event). 
For the CYGNUS array, the expected burst duration 
is no more than 40 ms for a source at zenith, while the 
dead time estimated from the data in Alexandreas et 
al. (1993) is <~ 30 ms per event. Since these two times 
are of the same order of magnitude, the real constraint 
on the PBH number density cannot be calculated. For 
the HEGRA array, we also disregarded the effect of its 
dead time, because no dead time is given in Funk et 
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al. (1995). Finally, estimates show that taking into 
account the dead time is unimportant at relatively low 
gamma-ray photon energies (~ 1 TeV). Therefore, the 
result of the Whipple array (Linton et al. 2006) re- 
mains essentially unchanged. 

Comparison of Figs. 6a and 6b leads us to conclude 
that taking into account the burst time profile and the 
array dead time is very important at high energies of 
the detected gamma-ray photons. 

How far are our constraints on the PBH number 
density from its actual estimates? If we make two as- 
sumptions: (1) the spectrum of the primordial density 
fluctuations is scale-invariant (in this case, the data ob- 
tained in measurements of the relic gamma radiation 
can be used to calculate the PBH mass spectrum and 
number density) and (2) the PBH accumulation factor 
in the Galaxy is equal to unity (i.e., there is no local 
increase in PBH number density in the Galaxy), then 
the predicted number density of evaporating PBHs is 
many orders of magnitude lower than the upper limits 
shown in Fig. 6b. However, there are many cosmo- 
logical models in which the spectrum of the primordial 
density fluctuations is nonmonotonic on small scales 
(see, e.g., Bugaev and Klimai (2006) and references 
therein). The possible value of the accumulation fac- 
tor has already been discussed above. In a word, we 
cannot rule out the possibility that the actual PBH 
number density is close to the already available exper- 
imental upper limits. In general, the predicted PBH 
number density in the Universe is determined by our 
views of the end of inflation. PBHs (even if they will 
not be discovered) are sources of unique information 
about this period in the history of the early Universe. 
Therefore, their searches will undoubtedly be contin- 
ued. 

As we see from Fig. 6b, the best (to date) con- 
straint on the number density of evaporating PBHs 
(1.08 x 10 6 pc -3 yr" 1 ) was obtained in the experi- 
ment on the Whipple Cherenkov telescope (Linton et 
al. 2006). However, it should be noted that the effec- 
tive gamma-ray photon energy in our experiment is two 
orders of magnitude higher than that in the Whipple 
one. Thus, our upper limit pertains not to black holes 
in general, but to black holes with certain properties 
(emitting 100- TeV gamma-ray photons at the end of 
their evaporation during bursts lasting ~ 10 ms). 
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